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The temporal evolution of the magnetization vector of a single-domain magnetostrictive nano- 
magnet, subjected to in-plane stress, is studied by solving the Landau-Lifshitz-Gilbert equation. 
The stress is ramped up linearly in time and the switching delay, which is the time it takes for the 
magnetization to flip, is computed as a function of the ramp rate. For high levels of stress, the 
delay exhibits a non-monotonic dependence on the ramp rate, indicating that there is an optimum 
ramp rate to achieve the shortest delay. For constant ramp rate, the delay initially decreases with 
increasing stress but then saturates showing that the trade-off between the delay and the stress (or 
the energy dissipated in switching) becomes less and less favorable with increasing stress. All of 
these features are due to a complex interplay between the in-plane and out-of-plane dynamics of the 
magnetization vector induced by stress. 



I. INTRODUCTION 

There is significant interest in studying the mag- 
netization reversal dynamics of multiferroic (strain- 
coupled magnetostrictive/piezoelectric bilayer) nano- 
magnets subjected to stress. This has potential applica- 
tions in ultralow power non-volatile magnetic logic and 
memory^— because switching a multiferroic nanomagnet 
with stress dissipates far less energy than switching it 
with a magnetic field or spin transfer torque produced 
by a currenti As a result, stress-mediated switching 
can reduce energy dissipation in magnetic reversal to the 
point where non-volatile memory and logic systems can 
be run by harvesting energy solely from the environment 
without needing a power source or battery! This opens 
up unique applications in situations where energy is a 
premium - such as in implanted medical devices, struc- 
tural health monitoring systems, "wearable" electronics, 
and space based applications. 

A multiferroic nanomagnet is made of a magnetostric- 
tive layer and a piezoelectric layer in intimate contact 
with each other (see Figure [1]). A voltage applied across 
the piezoelectric layer generates in it a mechanical strain 
that is mostly transferred to the magnetostrictive layer 
by elastic coupling and produces an extension^ if the 
latter layer is much thinner than the former. If we me- 
chanically constrain the magnetostrictive layer from ex- 
panding or contracting along a certain in-plane direction, 
e.g. along the minor axis of the ellipse in Figure [TJ then 
this will generate uniaxial stress along the major axis 
through the d^i coupling in the piezoelectric. This stress 
will cause the magnetization axis of the magnetostrictive 
layer (nanomagnet) to rotate by a large angle 10 , which 
has been demonstrated in recent experiments^, although 
not in nanoscale. 

Let us assume that the shape of the nanomagnet is that 




FIG. 1. An elliptical multiferroic nanomagnet stressed with 
an applied voltage. The polarity of the applied voltage V is 
determined by the sign of the magnetostrictive coefficient of 
the material used as the magnetostrictive layer. The polarity 
should be such that the stress generated favors aligning the 
magnetization vector along the in-plane hard axis rather than 
the easy axis. 



of an elliptical cylinder as shown in Figure [T] and that the 
initial orientation of the magnetization is close to the ma- 
jor axis of the ellipse (z-axis), which is the magnet's easy 
axis. In that case, a 90° rotation will place the magneti- 
zation vector along the minor axis, which is the in-plane 
hard axis. Subsequent removal of stress can relax it back 
to the easy axis, but in a direction anti-parallel to the 
initial direction, resulting in a ^180° rotation or "flip" 1 . 
Ref. |l| showed that the energy dissipated in this process 
is extremely small (~200 kT at room temperature) for 
optimum choice of materials, even when the switching 
takes place in ~1 nanosecond and the stress is turned on 
abruptly and instantaneously. In fact, a tiny voltage of 
few mV applied abruptly can generate enough stress to 
flip the magnetization^ in 1 nanosecond, which results 
in the miniscule dissipation of ^200 kT. 

In this paper, we are concerned with the following is- 
sue. The applied voltage cannot generate strain in the 
magnetostrictive layer instantaneously. If we ramp up 
the voltage gradually with a rise time longer than the 
response time of strain, then strain may be able to fol- 
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low the voltage quasi-statically. In that case, by con- 
trolling the ramp rate of the voltage, we can control the 
rise time of the strain. This may have significant effects 
on both the switching delay and the energy dissipated in 
the switching process. The purpose of this paper is to 
investigate this possibility. 

Intuitively, one would expect that if the stress is al- 
ways ramped up to a constant value regardless of the 
ramp rate, then the time taken to flip the magnetization 
(switching delay), will decrease monotonically with in- 
creasing ramp rate. The only caveat is that at high stress 
levels, a very fast ramp rate may cause ripples and ring- 
ing in the temporal evolution of the magnetization vector, 
which may prolong the switching process. The actual sit- 
uation turns out to be a little more complicated because 
the switching dynamics exhibits rich and complex be- 
havior as a result of the interplay between the in-plane 
and out-of-plane excursions of the magnetization vector 
under application of stress. This complex interplay has 
two effects: (1) it makes the switching delay exhibit a 
non-monotonic dependence on the ramp rate when high 
stresses are encountered, and (2) it makes the switching 
delay saturate quickly with increasing stress at a con- 
stant ramp rate. In this paper, we have studied this in- 
triguing dynamics by solving the Landau-Lifshitz-Gilbert 
(LLG) equation which governs the temporal evolution of 
the magnetization vector of a single domain nanomagnet. 

The rest of this paper is organized as follows. In Sec- 
tion [IIJ we first derive the torque exerted by any applied 
stress on the magnetization vector of the magnetostric- 
tive layer, and then solve the LLG equation analytically 
in the spherical coordinate system to yield equations that 
govern the time evolution of the magnetization vector. 
These equations describe how the polar angle 9{t) and the 
azimuthal angle cf>(t) of the magnetization vector change 
with time. They are solved numerically to obtain the dy- 
namics of magnetization rotation. In Section [ill! we de- 
rive the energy landscapes of the nanomagnet (energy as 
a function of magnetization orientation) in the stressed 
and relaxed conditions since they are a valuable aid in 
understanding the time evolution of the magnetization 
vector. In Section HVl we present the simulation results, 
while in Section [V] we discuss the implications of these 
results and present the conclusions. 



II. SOLUTION OF THE 
LANDAU-LIFSHITZ-GILBERT EQUATION 

Consider an isolated nanomagnet in the shape of an 
elliptical cylinder whose elliptical cross section lies in 
the y-z plane with its major axis aligned along the z- 
direction and minor axis along the y-direction. The di- 
mension of the major axis is a, that of the minor axis is 
b, and the thickness is The volume of the nanomag- 
net is 57 = (ir/A)abl. Let 6(t) be the angle subtended by 
the magnetization axis with the -t-z-axis at any instant 
of time t and 4>{t) be the angle between the +x-axis and 



the projection of the magnetization axis on the x-y plane. 
Thus, 8(t) is the polar angle and <fi(t) is the azimuthal 
angle. Note that when tfi = 90°, the magnetization vector 
lies in the plane of the magnet. Any deviation from (f> — 
90° corresponds to out-of-plane motion. 

The total energy of the single-domain nanomagnet 
(magnetostrictive layer) is the sum of the uniaxial shape 
anisotropy energy and the uniaxial stress anisotropy en- 
ergy: 



E{t) =E S HA(t)+E S TA(t), 



(1) 



where EsHA(t) is the uniaxial shape anisotropy energy 
and EsTA(t) is the uniaxial stress anisotropy energy at 
time t. The former is given by 



EsHAit) = {n /2)MfClN d {t), 



(2) 



where M a is the saturation magnetization and N d (t) is 
the demagnetization factor expressed asii 

N d (t) = N d - zz cos 2 0(t) + N d - yy sin 2 6{t) sin 2 <j>{t) 

+ N d 9(t) cos 2 (j)(t) (3) 

with N d — zz , N d ^ yy , and N d - xx are the components of the 
demagnetization factor along the z-axis, y-axis, and x- 
axis, respectively. When a and b are nearly equal, I a.b 



and a > b, N d _ zz , 
given bj*ii 



N d - yy , and N d - xx are approximately 
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More accurate expressions for these quantities can be 
found in Ref. 0. 

Note that uniaxial shape anisotropy will favor lining 
up the magnetization along the major axis (z-axis) by 
minimizing Esha, which is why we will call the major 
axis the "easy axis" and the minor axis (y-axis) the "hard 
axis" in the plane of the magnet. By mechanically con- 
straining the nanomagnet from expanding or contracting 
in the y-direction using appropriate clamps, we will gen- 
erate uniaxial stress along the z-axis (easy axis) . In that 
case, the stress anisotropy energy is given by 



E, 



STA 



(i) = -{3/2)X s a{t)ncos 2 9{t), 



(5) 



where (3/2)A s is the magnetostriction coefficient of the 
nanomagnet and a(t) is the stress generated in it by an 
external agent. Note that a positive A s <r(t) product will 
favor alignment of the magnetization along the major 
axis (z-axis), while a negative X s a(t) product will favor 
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alignment along the minor axis (y-axis), because that 
will minimize Est A (t) ■ In our convention, a compressive 
stress is negative and tensile stress is positive. Therefore, 
in a material like Terfenol-D that has positive \ s , a com- 
pressive stress will favor alignment along the minor axis, 
and tensile along the major axis. The situation will be 
opposite with nickel and coblat that have negative A s . 

At any instant of time, the total energy of the nano- 
magnet can be expressed as 

E(t) = E{6{t), <j>{t)) = B(t)sin 2 6(t) + C(t) (6) 

where 

B (t) = B (ct>{t)) = yM s 2 fi [N d _ xx cos 2 4(t) 

+N d -yysin 2 <fi(t) — Nd-zz] (7a) 

B st ress(t) = (3/2)A s CT(i)fi (7b) 

B(t) = B (t) + B 

stress 

(t) (7c) 
C(t) = ^M 2 ilN d ^zz - (3/2)A.<r(t)fi. (7d) 

Note that Bo(t) is always positive for our choice of ge- 
ometry, but B stre ss{t) can be negative or positive in ac- 
cordance with the sign of the \ s a(t) product. 

The magnetization M(t) of the nanomagnet has a con- 
stant magnitude at any given temperature but a variable 
direction, so that we can represent it by the vector of unit 
norm n m (£) = M(t)/|M| = e r where e r is the unit vec- 
tor in the radial direction in spherical coordinate system 
represented by (r,9,(f>). The other two unit vectors in the 
spherical coordinate system are denoted by eg and for 
9 and (/) rotations, respectively. The gradient of potential 
energy at any particular instant of time t is given by 

VE(t) = VE(6{t),ct>{t)) = ^Sr&e + —X^^Sle, 



dO(t) sinO(t) d<t>{t) 



where 

dE{t) 
d6{t) 



(8) 



2B(t)sin9{t)cos9(t) - (3/2)\ s Qcos 2 9(t) 



2 x 5<7 (*) 
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= 2B(t)sin0(t)cos8(t) - cr c (t) 



da(t) 

mt) 



(9) 



9E(t) _ 
d<j){t) 2 



M 2 n(N d _ xx - N d _ yy )sin(2(b(t)) S in 2 9(t) 



-(3/2)\ s ncos 2 6(t) d<T{ " 



= -B 0e (t) sin 2 9{t) - <r B (t) 



dm 

da(t) 

dW) 



(10) 



while 



B 0e (t) = B 0e m)) 



M 2 n{N d . 



a c (t) = (3/2)X s ncos 2 d(t). 



N d - yy )sin{2(t>(t)) 
(11) 



The torque acting on the magnetization per unit volume 
due to shape and stress anisotropy is 



T B (t) = -n m (t) x V£(0(i), <£(*)) 
= e r x 



2B(t)sin6(t)cos0(t) - a c (t)^±) e e 
2B(t)sin6(t)cos9(t) - °c(t)^^j e^ 

Boi(() , i „ e(t) + _ S |_||j s „ (12) 

This torque causes the magnetization vector to rotate. 
The magnetization dynamics under the action of this 
torque is described by the Landau-Lifshitz-Gilbert (LLG) 
equation as follows. 



dn m (t) ( dn m (t) 
a n m (t) x 



df 



dt 



7 
M v 



■T B (t) (13) 



where a is the dimensionless phenomenological Gilbert 
damping constant, 7 = 2/iB/i /?i is the gyromag- 
netic ratio for electrons and is equal to 2.21 x 10 5 
(rad.m).(A.s)" 1 , \ib is the Bohr magneton, and My — 
HqM s VL. In the spherical coordinate system, 



dn m (t) dd(t) 



dt 



dt 



sin9{t) — 



(14) 



Accordingly, 



( '< ( o m (7) X dn ™^ j = -nsintiij )</(/) do + (\ti'{t)e<f > 

(15) 

where ()' denotes d()/dt. This allows us to write 



dt 



dn m (t) ( dn m (t) 
a n m (t) x 



dt V dt 

= (9'(t)-asin9(t) 4>'{t)) e + {sin9(t) 4>'(t)+a9'(t)) e . 

(16) 

Equating the eg and components in both sides of 
Equation ((T3|) . we get 



-asin9(t) (j>'{t) 

_ J_ 

M v 



sin0(t) (f>'(t)+a6'(t) = 

- jt- (2B(t)sind{t)cos6(t) - *c(oS|) . (18) 

Solving the above equations, we get the following coupled 
equations for the dynamics of 9(t) and cf>(t). 



4 



d9{t) 
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<l + a 2 ) 



dt 
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dt 
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(19) 
(20) 



We should note that the Equations [[T2" ]) . ([19 ]) . and ^ 
are not valid when sin0 = (0 = 0° or = 180°), 
i.e. when the magnetization direction is exactly along 
the easy axis. At these points, the torque on the mag- 
netization vector given by Equation f|12[) diverges. In 
order to avoid these points, we will assume that the ini- 
tial orientation of the magnet is 6 — 179°, and switching 
is deemed to have been completed when = 1°. This 
1° deflection could be caused by thermal fluctuations. 
Similar assumptions have been made by other authors^. 

Equation (| 12|) shows that there is an internal feedback 
in the system. Stress induces a torque which produces 
the rotation (0(f), <i>{t)). That rotation generates an addi- 
tional torque through the da/88 and da /dip terms. That 
additional torque affects the response. This feedback 
mechanism determines the relation between the rotation 
and stress, and hence the switching delay as a function 
of stress. 

Note from Equation (TT2"j) that the torque has contribu- 
tions due to the dynamic change in stress (da(t)/d8(t), 
da(t)/d(j)(t) terms). These contributions may aid or hin- 
der the rotation of the magnetization vector at differ- 
ent times. This is why the switching delay will depend 
on the ramp rate da/dt. This dependence turns out to 
be non-monotonic because of the complex actions of the 
magnetization vector. 

Note from Equation (TTTj) that the term Bo e (t) will be 
negative when 90° < </> < 180°. In that case, its con- 
tribution to the time rate of change of 0, i.e. d0(t)/dt, 
will be positive, as we can see from Equation (fT9"|) . In 
other words, a negative Bo e (t) will tend to increase 
with time. Since our initial value of is 179° and the final 
value is 1°, we would prefer that will always decrease 
- and never increase - with time in order to complete 
the switching in the shortest time. Therefore, a negative 
value of Bo e (t), or equivalently <f> lying in the interval 
[90°, 180°], is counterproductive since that makes in- 
crease with time, causing the magnetization to rotate in 
the wrong direction, opposite to the preferred direction. 
This hinders switching and increases the switching delay. 
Therefore, we will always prefer that Bo e (t) remains pos- 
itive, or equivalently remains in the interval [0°, 90°] 
or [180°, 270°]. 

On the other hand, it is clear from Equation (|2T)|) that 
a positive Bo e (t) makes a positive contribution to the 
rate dcj>/dt, which will tend to increase <fi with time and 
make it exceed 90°. These two counteracting influences 
of BQ e (t) determine the actual switching dynamics and 
the resulting switching delay. 
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FIG. 2. Energy landscape of a Terfenol-D/PZT multiferroic 
nanomagnet. Plot of total energy as a function of polar and 
azimuthal angles of the magnetization vector. 

Another point to note is that when the applied stress 
is sufficiently high, the stress term B stre ss(t) dominates 
the term B(t) in Equations (TJSJ and ((101). The term in- 
volving B{t) should remain positivein Equation (|19p in 
order to help dd{t) / dt remain negative so that the magne- 
tization vector can rotate in the right direction. In order 
to keep B stre ss{t) negative, we will have to ensure that 
the product X s a is negative. For materials with positive 
magnetostriction (e.g. Terfenol-D), this requires that the 
stress is negative, while for materials with negative mag- 
netostriction (e.g. nickel or cobalt), the stress should 
be positive. Since tensile stress is positive and compres- 
sive is negative, Terfenol-D will require compressive stress 
and nickel or cobalt will require tensile stress to initiate 
switching if the magnetization is initially aligned close to 
the easy axis. 

III. ENERGY LANDSCAPE 

The energy landscape of a nanomagnet, which plots 
the total energy E(t) as a function of the polar angle 0(f) 
and the azimuthal angle (f>(t), provides valuable informa- 
tion. The final state of the magnetization will always be 
at an energy minimum. Stress will modify the energy 
landscape of a nanomagnet and shift the energy mini- 
mum from one set of angles (4>i,9i) to another (0/,0/), 
thereby effecting switching of the magnetization. 

Figure [5] shows the energy landscape of a Terfenol- 
D/PZT multiferroic nanomagnet for 0° < < 180° and 
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Energy landscape with +=0 degree (TerfenolD) 




e (degree) 
00 

FIG. 3. Energy landscape (E versus 9) of a Terfenol-D/PZT 
multiferroic nanomagnet for (a) = 0° and (b) <j> = 90° . The 
nanomagnet parameters are given in Table I. The maximum 
stress that can be generated in a Terfenol-D/PZT multifer- 
roic nanomagnet is of the order of 40 MPa. This stress is 
enough to depress the shape anisotropy barrier and switch 
the magnetization from the easy axis (8 = 0°, or 180°) to the 
in-plane hard axis (9 = 90°) by shifting the energy minimum 
to 9 — 90° when the magnetization vector lies in the plane 
of the magnet. The shape anisotropy energy barrier is much 
taller when cj> = 0° than when cj> = 90° because of the small 
thickness of the nanomagnet making Nd- xx much larger than 

N d -yy Or Nd-ZZ- 



in the magnetostrictive layer by the strained piezoelec- 
tric layer is usually sufficient to rotate the magnetization 
axis when <j) = 90°, because the shape anisotropy bar- 
rier that has to be overcome is small. However, this does 
not necessarily happen when <fi — 0° because then the 
shape anisotropy energy barrier is much higher. On the 
other hand, out-of-plane excursion of the magnetization 
vector generates an additional torque that aids switch- 
ing. Thus, some out-of-plane excursion is beneficial, but 
if the magnetization vector strays out of plane, it en- 
counters a larger energy barrier that prevents switching. 
Therefore, the magnetization vector must precess and ul- 
timately return close to the nanomagnet's plane before 
switching can be accomplished. This is the cause of pre- 
cessional motion. 

The energy landscapes allow us to estimate the min- 
imum stress needed to rotate the magnetization vector 
from the initial orientation close to the easy axis to the 
in-plane hard axis. Once the vector aligns along the in- 
plane hard axis, stress is removed. Thereafter, the mag- 
netization vector will relax back to the easy axis but to an 
orientation anti-parallel to the initial orientation. This 
results in switching. The minimum stress required for 
this purpose is found by equating the shape anisotropy 
energy barrier to the stress anisotropy energy, i.e. 

(N d - yy - N d _ zz ) = ^X s a, (21) 

which yields 

Omin = {N d -yy - N d - zz ) M 2 S ^-. (22) 

However, switching with this minimum stress will incur a 
very long switching delay, so that some excess stress will 
be needed to switch the magnetization reasonably fast. 
Normally, one would expect the switching delay to de- 
crease continuously with increasing excess stress, but in 
reality it saturates beyond a certain stress so that increas- 
ing stress further offers only marginal advantage. This 
feature cannot be understood from the energy landscape 
because it is a consequence of the interplay between the 
in-plane and out-of-plane dynamics of the magnetization 
vector, which is not captured in the energy profiles. 



0° < <\> < 360° without any applied stress. Figures 3(a) 
and |3(b)| show how stress modifies the energy landscape 
in 0-space for = 0° and = 90°, respectively. 

When = 0°, i.e. the magnetization vector lies in 
the x-z plane, the energy barrier separating the two sta- 
ble magnetizations along the z-axis (easy-axis) is ~10 
times taller than what it is when <f> = 90°, i.e. when the 
magnetization vector lies in the y-z plane. This happens 
because the thickness of the nanomagnet is much smaller 
than the other two dimensions, which makes the shape 
anisotropy energy barrier much taller in the former case 
than in the latter case. The stress that can be generated 



IV. SIMULATION RESULTS 

We consider a multiferroic nanomagnet composed of a 
PZT layer (lead-zirconate-titanate) and a magnetostric- 
tive layer which is made of either polycrystalline Terfenol- 
D, or polycrystalline nickel, or polycrystalline cobalt. Be- 
cause it is polycrystalline, the magnetocrystalline layer 
does not have significant magnetocrystalline anisotropy. 
The material parameters for the magnetostrictive layer 
are given in Table HJ 4 ~ 1S ■ They ensure that the shape 
anisotropy energy barrier is ^32 kT at room tempera- 
ture. The PZT layer is assumed to be four times thicker 
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TABLE I. Material parameters for different materials. 





Terfenol-D 


Nickel 


Cobalt 


Major axis (a) 


101.75 nm 


105 nm 


101.75 nm 


Minor axis (b) 


98.25 nm 


95 nm 


98.25 nm 


Thickness (t) 


10 nm 


10 nm 


10 nm 


Young's modulus (Y) 


8xl0 lu Pa 


2.14xlO ii Pa 


2.09X10 11 Pa 


Magnetostrictive coefficient ((3/2)A s ) 


+90xl0~ s 


-3xl0~ b 


-3xl0~ s 


Saturation magnetization (M s ) 


8xlO b A/m 


4.84xlO b A/m 


8xlO b A/m 


Gilbert's damping constant (a) 


0.1 


0.045 


0.01 



R 




V 



t=0 t=T t=T end 

FIG. 4. The switching circuit and ramp parameters. 



time considered in our simulation is 1-150 ps, we are in 
the adiabatic limit (T ^> RC) and hence the energy dis- 
sipation in the external circuit that generates the voltage 
V across the PZT layer will be less than (1/2)CV 2 . We 
assume that the charging circuit is represented by the cir- 
cuit diagram in Figure |4] The energy dissipated -Ed, rise 
during the rise of the voltage (charging cycle) for a sig- 
nal of total time-period T en d and ramp-period T can be 
calculated as below. 



than the magnetostrictive layer so that any strain gen- 
erated in it is transferred almost completely to the mag- 
netostrictive layer. We will assume that the maximum 
strain that can be generated in the PZT layer is 500 
ppmi^, which would require a voltage of 1 1 1 mV because 
<i3i=1.8e-10 m/V for PZT—. The corresponding stress is 
the product of the generated strain (500 x 10~ 6 ) and the 
Young's modulus of the magnetostrictive layer. Based 
on available data for Young's modulus, the maximum al- 
lowable stresses for Terfenol-D, nickel, and cobalt are 40 
MPa, 107 MPa, and 104.5 MPa, respectively. 

In all our simulations, the initial orientation of the 
magnetization vector is: 9 = 179° and (j> — 90°. Stress 
is applied as a linear ramp and we solve Equations (|19[) 
and (|20|) at each time step. Once 9 becomes 90°, stress is 
removed and we follow the magnetization vector in time 
until 9 becomes 1°. At that point, switching is deemed 
to have occurred. 

We assume that the voltage applied on the piezoelectric 
is ramped up linearly to its steady-state value in time T 
which we call the rise time. When the stress is ramped 
down, we use the same rate, i.e. we reduce the stress 
from its maximum value to zero in time T . In all cases, 
the rise time is equal to the fall time. 

We also assume that the PZT layer, which acts as a 
capacitor, is electrically accessed with a silver wire of re- 
sistivity ^2.6 /ifi-cm^i so that a typical access line of 
length 10 nm and cross section 50 nm x 50 nm will have 
a resistance of ~100 f2. Based on the dimensions of the 
PZT layer (major axis, minor axis, and thickness), and 
assuming that the relative dielectric constant of PZT is 
1000, the capacitance of the PZT layer will be -2 fF. 
Therefore, the RC time constant associated with charg- 
ing the capacitor is ~0.2 ps. Since the range of ramp 



E, 



d,rise 



/ RC\ ( RC RC rr i D/ ~. 

= cv 2 — 1 - — + — e - T/RC 



1 (RC 
~2 V~ 



V T J I T ' T 

1 _ e -T/RCy e -2(T end -T)/RC y 



where C is the capacitance of the PZT layer and V is the 
steady state voltage that generates the required stress. 
The last term in the above expression comes from a finite 
value of T end . 

The energy dissipated during the discharging cycle is 
Edjaiii which can be calculated from an expression simi- 
lar to the one above, except that the value of T end may be 
different. For the sake of brevity, we will term the total 
energy dissipated in the charging circuit E d , r ise + E d j a u 
as the l CV 2 ' energy dissipation in the remainder of this 
paper. 

Because of Gilbert damping in the magnet, an ad- 
ditional energy Ed is dissipated when the nanomagnet 
switches. This energy is given by the expression 



E,, 



Pd(t)dt, 



(24) 



where r is the switching delay and Pd(t), which is the 
dissipated power, is given b y 22 ' 23 



Pd(t) 



a 7 



(1 + a 2 )/i Af s fi 



mt)\ 



(25) 



We sum up the power Pd{t) numerically throughout 
the switching period to get the corresponding energy dis- 
sipation Ed and add that to the 'CV 2 ' dissipation in 
the switching circuit to find the total dissipation E tota \. 
The average power dissipated during switching is simply 

E d /T. 
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Dynamics of 8 with time (40 MPa stress, 1 ps ramp) 
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FIG. 5. Magnetization dynamics in the Terfenol-D/PZT mul- 
tiferroic nanomagnet. The stress is ramped up from to 40 
MPa in 1 ps: (a) polar angle 9 versus time, and (b) the tra- 
jectory traced out by the tip of the magnetization vector in 
three-dimensional space while switching occurs, i.e. during 
the time 9 changes from 179° to 1°. Note that the magnet's 
plane is x = 0. 

We analyze the magnetization dynamics of the magne- 
tostrictive layer as a function of both the magnitude and 
the ramp time of the stress for three different materials 
(Terfenol-D, nickel, cobalt). The results are presented 
in the ensuing subsections. These materials are chosen 
because of their very different material parameters such 
as Gilbert damping constant, saturation magnetization, 
and magnetostrictive coefficient. 

A. Terfenol-D 

Terfenol-D has a positive magnetostrictive coefficient 
(see Table U). Therefore, we will need a compressive 
stress to rotate the magnetization vector away from its 
initial alignment close to easy axis (6 = 179°) towards 
9 = 90°. Note that we need to use the correct voltage 
polarity to ensure that a compressive stress is generated 
on the Terfenol-D layer. The maximum stress that can 
be generated on the Terfenol-D layer with the maximum 
allowed 500 ppm strain in the PZT layer is 40 MPa, and 



Dynamics of e with time (40 MPa stress, 60 ps ramp) 
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FIG. 6. Magnetization dynamics in the Terfenol-D/PZT mul- 
tifcrroic nanomagnet. The stress is ramped up from to 40 
MPa in 60 ps: (a) polar angle 9 versus time, and (b) the 
trajectory traced out by the tip of the magnetization vector 
while switching occurs. 

the minimum stress that is needed to switch the nano- 
magnet is found by equating the stress anisotropy energy 
to the shape anisotropy energy barrier. This stress is 1.91 
MPa. 



1. Ramp rate and switching delay 

Equations (fl"9| and ([20]). derived in Section lU are 
solved numerically to find the values of 9{t) and </>(£) 
at any given instant t. This yields the magnetization 
dynamics under various stresses and ramp rates. 

a. Fast ramp: The stress on the Terfenol-D layer is 
ramped up linearly in time from to the maximum pos- 
sible value of 40 MPa in 1 ps. The corresponding mag- 
netization dynamics is shown in Figure [5] We notice 
that the polar angle continuously evolves from its ini- 
tial value of 179° towards its final value of 1° for the first 
200 ps. However, because of the coupled dynamics of the 
azimuthal angle cf>(t) and the polar angle 6(i), 4>(t) devi- 
ates from its initial value of 90° at around 200 ps, forcing 
the magnetization vector to venture out of plane. This 
makes the magnetization vector execute precessional mo- 
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Ramp time vs. switching delay (Terfenol-D) Stress vs. switching delay (Terfenol-D) 




Ramp time (ps) Stress (MPa) 



FIG. 7. Switching delay of the Terfenol-D/PZT nanomagnet 
as a function of the rise (or fall) time of the ramp, with the 
magnitude of stress as a parameter. 



FIG. 10. Switching delay versus stress for the Terfenol- 
D/PZT multiferroic nanomagnet for two different ramp rise 
(and fall) times of 1 ps and 150 ps. 



Dynamics of e with time (40 MPa stress, 40 ps ramp) 
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FIG. 8. Magnetization dynamics in the Terfenol-D/PZT mul- 
tiferroic nanomagnet when the stress is ramped up linearly 
from to 40 MPa in 40 ps. The polar angle is plotted ver- 
sus time. There is a slight ripple, but its amplitude is greatly 
reduced compared to the case when the rise (and fall) time is 
1 ps. The switching delay in this case is about 285 ps. 
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FIG. 9. Magnetization dynamics in the Terfenol-D/PZT mul- 
tiferroic nanomagnet when the stress is ramped up linearly 
from to 40 MPa in 50 ps. The polar angle 9 is plotted ver- 
sus time. There is again a slight ripple, but its amplitude is 
very small and barely perceptible. The switching delay in this 
case is about 500 ps. 



FIG. 11. Energy dissipated in flipping the magnetization of 
the Terfenol-D/PZT multiferroic nanomagnet as a function 
of switching delay for a ramp rise (and fall) time of 150 ps. 
This range of switching delay corresponds to a stress range of 
2.5 MPa to 40 MPa. The energy dissipated in the nanomag- 
net due to Gilbert damping and the energy dissipated in the 
external switching circuit ('CV 2 ') are shown separately. 
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FIG. 12. For a fixed stress of 40 MPa, energy dissipated in 
flipping the magnetization of the Terfenol-D/PZT nanomag- 
net as a function of switching delay when the latter is varied 
by varying the ramp's rise (and fall) time from 80-150 ps. 
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Stress: 10 MPa, Ramp Time: 1-150 ps (Terfenol-D) 
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FIG. 13. For a fixed stress of 10 MPa, energy dissipated in 
flipping the magnetization of the Terfenol-D/PZT nanomag- 
net as a function of switching delay when the latter is varied 
by varying the ramp's rise (and fall) time from 1-150 ps. 



Stress: 1.91 MPa, Ramp Time: 1-150 ps (Terfenol-D) 
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FIG. 14. For a fixed stress of 1.91 MPa, energy dissipated in 
flipping the magnetization of the Terfenol-D/PZT nanomag- 
net as a function of switching delay when the latter is varied 
by varying the ramp's rise (and fall) time from 1-150 ps. 



tion in space while its projection on the magnet's plane 
changes course and rotates in the direction opposite to 
the desired direction so that 9(t) begins to increase with 
time instead of decreasing. Eventually, as the magnetiza- 
tion vector stops precessing and returns to the magnet's 
plane, which happens around 330 ps, its projection on 
the magnet's plane starts rotating towards its final desti- 
nation, ultimately reaching 9 = 1° and <\> — 90°. Because 
of the interplay between the 9- and ^-dynamics, which 
causes the magnetization vector to leave the plane of the 
nanomagnet after ^200 ps, switching takes around 700 
ps. 

b. Slow ramp: Figure [6] shows the magnetization dy- 
namics for a slow ramp that takes 60 ps to rise linearly 
from to 40 MPa. In this case, the magnetization vec- 
tor does not even budge from its initial orientation of 
9 = 179° until the stress reaches its peak value of 40 
MPa, which happens after 60 ps have elapsed. There- 
after, the magnetization vector rotates towards its final 
destination of 9 = 1° without ever changing course and 



rotating in the opposite direction, unlike the previous 
case. The magnetization vector however does leave the 
magnet's plane in this case as well, but clearly it does 
not precess as much as in the previous case (see the tra- 
jectory plot). The switching is actually faster now and 
takes only 250 ps compared to 700 ps for the previous 
case. Thus, a slower ramp can be beneficial when high 
stresses are applied. It eliminates the ripple and ringing 
in the switching characteristic seen in the previous case 
by limiting the out-of-plane excursion of the magnetiza- 
tion vector and saves precious time. 

Figure [7] shows the switching delay as a function of the 
ramp's rise (or fall) time with the magnitude of stress 
as a parameter. We see that for stresses of 5 MPa and 
10 MPa, the switching delay increases linearly with the 
ramp time, but for higher stresses of 30 MPa and 40 
MPa, the switching delay shows clear non-monotonic be- 
havior for rise (and fall) times less than 60 ps. Normally, 
the switching delay should increase continuously with the 
ramp's rise (and fall) time, but the out-of-plane dynam- 
ics and precession of the magnetization vector that occurs 
at very fast ramp rates can reverse this trend and cause 
the non-monotonic behavior. Figures [5] and [5] show the 
magnetization dynamics for ramp rise (and fall) times of 
40 ps and 50 ps, respectively. There are some ripples in 
these cases also, but their amplitudes are much reduced 
so that they are barely visible in the plots. For ramp rise 
times exceeding 60 ps, the ripples disappear. Thereafter, 
the switching delay increases monotonically with the rise 
(or fall) time. 



2. Switching delay and energy dissipation 

Figure [TO] shows the dependence of the switching delay 
on stress in the Terfenol-D/PZT nanomagnet, with the 
ramp's rise (or fall) time as a parameter. The two rise 
(and fall) times considered are 1 ps and 150 ps. There 
is a cross over at 14 MPa stress. At low stress levels 
below 14 MPa, not much ripple is generated by a fast 
ramp so that the switching delay is shorter for the faster 
ramp. At high stress levels exceeding 14 MPa, a fast 
ramp generates enough ripple that the switching delay 
becomes longer for the faster ramp. This is the reason 
for the cross over. 

Figure [IT] shows the energy dissipated in flipping the 
magnetization of the Terfenol-D/PZT multiferroic nano- 
magnet as a function of switching delay for a fixed rise 
(and fall) time of 150 ps. The switching delay is varied 
by varying the stress on the nanomagnet between 2.5 and 
40 MPa. The energy dissipated internally in the nano- 
magnet (Ed) and the energy dissipated in the switching 
circuit ('CV 2 ') are shown separately. They both tend to 
saturate at larger delays. 

For longer switching delays, the stress needed to flip 
the magnetization is less and hence the voltage V needed 
to generate the stress is smaller. This leads to a smaller 
'CV 2 ' dissipation in the switching circuit. At the same 
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time, the energy Ed dissipated internally in the nano- 
magnet is smaller when we switch slowly. In this range 
of switching delay, the energy dissipated in the external 
circuit is much smaller than the energy dissipated inter- 
nally in the nanomagnet since the switching is adiabatic 
(the rise (and fall) time is much longer than the RC time 
constant of the switching circuit). The ratios of the two 
energies however decreases with decreasing switching de- 
lay. Below a switching delay of 4 ns, the energy dissipated 
internally in the nanomagnet and energy dissipated in the 
switching circuit both increase super-exponentially with 
decreasing switching delay. At 1 ns switching delay, the 
total energy dissipated in switching is only about 200 
kT, which makes this switching methodology extremely 
energy- efficient. For this switching delay, the energy dis- 
sipated to switch a state-of-the-art transistor would have 
been at least two orders of magnitude larger—, and the 
energy dissipated to switch the same nanomagnet with 
spin transfer torque will be also at least two orders of 
magnitude larger—. 

Figure [12] plots the energy dissipation as a function of 
switching delay for a fixed stress of 40 MPa. Here, the 
switching delay is varied by varying the ramp's rise (and 
fall) time between 80 and 150 ps. A lower rise (or fall) 
time is avoided because of the non-monotonic behavior 
observed in Figure[7J Both Ed and the l CV 2 ' dissipations 
fall off with increasing switching delay. This means that 
the average power dissipation falls off even more rapidly 
with increasing switching delay since the energy dissipa- 
tion is the product of the average power dissipation and 
the switching delay. This figure shows that we can switch 
in 0.3 - 0.4 ns by dissipating 600 - 700 kT of energy at 
room temperature. 

A similar plot for a lower fixed stress of 10 MPa is 
shown in Figure 1131 Here, the switching delay is var- 
ied by varying the rise (and fall) time between 1 ps and 
150 ps since there is no issue of non-monotonic behav- 
ior at such low stress values (no ripples generated in the 
switching characteristics). We see that the energy dissi- 
pated internally in the nanomagnet (Ed) now increases 
with increasing switching delay which is the opposite of 
the behavior observed in the case of 40 MPa stress. In 
this case, the average power dissipation still goes down 
with increasing switching delay, but not fast enough, so 
that the total energy, which is the product of the average 
power and the switching delay, actually go up with in- 
creasing delay. However, the l CV 2 ' energy dissipation in 
the external circuit does goes down with increasing de- 
lay since switching becomes more adiabatic as the delay 
becomes longer. This figure shows that we can switch in 
~1 ns by dissipating roughly 200 kT of energy. 

Figure [H] shows the energy dissipation as a function 
of switching delay for a fixed stress of 1.91 MPa. Once 
again, the delay is varied by varying the ramp's rise (and 
fall) time between 1 and 150 ps. In this case, Ed is 
nearly independent of the switching delay, meaning that 
the average power dissipation in the nanomagnet varies 
inversely with the switching delay. The 'CF 2 ' dissipation 



in the external circuit still goes down with increasing de- 
lay as expected because switching becomes increasingly 
'adiabatic'. This figure shows that we can switch in 110 
ns by dissipating only ^35 kT of energy at room temper- 
ature. The corresponding average power dissipation in 
this case is roughly 35 kT/110 ns = 1.33 pW per nano- 
magnet per bit flip. If we have an array of magnets with 
areal density 10 10 cm~ 2 (10 Gbits/sq-cm) and 10% of 
them are being flipped at any given time (10% activity 
level), then the power dissipated is 1.3 mW/cm 2 . The 
energy needed to run at such low power levels can be har- 
vested from the local surroundings using existing energy 
harvesting devices 2 ^— without requiring a battery or ex- 
ternal energy source! This opens up the possibility of 
unique applications such as medically implanted devices 
(e.g. processors implanted in a patient's brain which 
warn of impending epileptic seizures) that run by har- 
vesting energy from a patient's body movements without 
requiring a battery, buoy mounted processors in the open 
sea that harvest energy from swaying motion induced 
by sea waves, or distributed sensor-processor networks 
for structural health monitoring of bridges and build- 
ings that harvest energy from vibrations of the structure 
due to wind or passing traffic. These applications are 
made possible by the extreme energy efficiency of strain- 
induced switching. 



B. Nickel 

Nickel has a negative magnetostrictive coefficient (see 
Table [TJ) so that a nickel/PZT multiferroic nanomagnet 
will require a tensile stress to initiate rotation away from 
the easy axis. However, the torque generated due to the 
change in stress with time (see Equations (IT91 and (|2T)]) ) 
would be of same sign as that in the case of Terfenol-D 
since it depends on the product of the magnetostrictive 
coefficient (A s ) and the change in stress (da). For the di- 
mensions of the nanomagnet chosen, the minimum stress 
that we will need in a nickel/PZT multiferroic to switch 
is 57 MPa, while the maximum stress that can be gen- 
erated by the 500 ppm strain in the PZT layer is 107 
MPa. 



1. Ramp rate and switching delay 

Just as in the case of Terfenol-D, Equations p^|) 
and (|20[) are solved numerically to find the values of 0(t) 
and 4>(t) at any given instant t for the nickel/PZT nano- 
magnet. This yields the magnetization dyanmics under 
various stresses and ramp rates. 

a. Fast ramp: Figure [15] shows the magnetization 
dynamics of a nickel/PZT multiferroic nanomagnet when 
the stress is ramped up linearly in time from to the 
maximum possible value of 107 MPa in 1 ps. The in- 
plane and out-of-plane dynamics of the magnetization 
vector are very similar to the case of Terfenol-D, except 
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Dynamics of e with time (107 MPa stress, 1 ps ramp) 
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FIG. 15. Magnetization dynamics in the nickel/PZT multi- 
ferroic nanomagnet. The stress is ramped up from to 107 
MPa in 1 ps: (a) polar angle 8 versus time, and (b) the tra- 
jectory traced out by the tip of the magnetization vector in 
three-dimensional space while switching occurs, i.e. during 
the time 9 changes from 179° to 1°. 

now we see a more ripples since there is more out-of- 
plane precession as can be seen in Figure [TST b). Nickel 
shows more ripples and more precession because it has a 
smaller Gilbert clamping constant than Terfenol-D. Con- 
sequently, the precessional motion is less damped. 

b. Slow ramp: Figure [TBT b) shows that a slow ramp 
rate (150 ps rise (and fall) time) quenches the preces- 
sion of the magnetization vector. This happens because 
a slower ramp causes less out-of-plane excursion of the 
magnetization vector and hence less precessional motion. 
Note that contrary to expectation, the slow ramp (150 ps 
rise (and fall) time) completes the switching in 0.95 ns, 
while the fast ramp (1 ps rise (and fall) time) takes twice 
the amount of time, i.e. 1.9 ns! This happens because 
of the ripples and ringing which have a deleterious effect 
on switching speed. Note that both ramp rates make 6 
come very close to 1° in about 0.9 ns, but the faster ramp 
rate causes the magnetization vector to pull back from 
the final destination and vacillate before reaching the fi- 
nal destination. This causes the ringing which prolongs 
the switching duration and increases the delay. 

Figure [T7] shows switching delay as a function of the 



Dynamics of e with time (107 MPa stress, 150 ps ramp) 
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FIG. 16. Magnetization dynamics in the nickel/PZT multi- 
ferroic nanomagnet. The stress is ramped up from to 107 
MPa in 150 ps: (a) polar angle 6 versus time, and (b) the 
trajectory traced out by the tip of the magnetization vector 
while switching occurs. 

ramp's rise (and fall) time for various stresses. At higher 
stress levels (80 - 107 MPa), the switching delay actually 
decreases with increasing rise (and fall) time (a counter- 
intuitive result) because of the ripples that are generated 
when the stress is ramped up to high values very fast. 
Thus, the out-of-plane dynamics of the magnetization 
vector plays a vital role in switching. 

For lower stress levels (60 -70 MPa), the switching de- 
lay increases slightly with increasing rise (or fall) time. 
Low stresses do not cause significant out-of-plane excur- 
sion of the magnetization vector leading to precessional 
motion. Therefore, even a fast ramp does not cause too 
many ripples and delay the switching when the stress is 
sufficiently low. In that case, a slower ramp will lead to 
slightly slower switching and that is what we observe. 

Strangely, at slow ramp rates (e.g. rise (and fall) time 
= 150 ps) and low stress levels, a slightly lower stress can 
cause considerably more precession than a slightly higher 
stress which leads to very strong dependence of switching 
delay on the magnitude of stress at low stress levels. In 
order to illustrate this, we plot the switching dynamics 
for 150 ps rise (and fall) time in Figures [18] and [19] for 
stresses of 60 MPa and 80 MPa. There are ripples and 
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Ramp time vs. switching delay (nickel) 
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FIG. 17. Switching delay of the nickel/PZT nanomagnet as 
a function of the rise (or fall) time of the ramp, with the 
magnitude of stress as a parameter. 
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FIG. 18. Magnetization dynamics in the nickel/PZT multifer- 
roic nanomagnet when the stress is ramped up linearly from 
to 60 MPa in 150 ps - (a) The polar angle 8 is plotted ver- 
sus time. There is some ripple, but its amplitude is reduced 
compared to the case when the rise (and fall) time is 1 ps. 
The switching delay in this case is about 6 ns. (b) trajectory 
traced out by the tip of the magnetization vector while the 
switching occurs. 
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FIG. 19. Magnetization dynamics in the nickel/PZT multifer- 
roic nanomagnet when the stress is ramped up linearly from 
to 60 MPa in 150 ps - (a) The polar angle 9 is plotted versus 
time. There is some ripple, but its amplitude is reduced com- 
pared to the case when the rise (and fall) time is 1 ps. The 
switching delay in this case is about 1.6 ns. (b) trajectory 
traced out by the tip of the magnetization vector while the 
switching occurs. 
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Stress vs. switching delay (nickel) 
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FIG. 20. Switching delay versus stress for the nickel/PZT 
multiferroic nanomagnet for two different ramp rise (and fall) 
times of 1 ps and 150 ps. 
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Stress: 60-107 MPa, Ramp time: 150 ps (nickel) 
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FIG. 21. Energy dissipated in flipping the magnetization 
of the nickel/PZT multiferroic nanomagnet as a function of 
switching delay for a ramp rise (and fall) time of 150 ps. This 
range of switching delay corresponds to a stress range of 60 
MPa to 107 MPa. The energy dissipated in the nanomag- 
net due to Gilbert damping and the energy dissipated in the 
external switching circuit ('CV 2 ') are shown separately. 
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FIG. 22. For a fixed stress of 107 MPa, energy dissipated 
in flipping the magnetization of the nickel/PZT nanomagnet 
as a function of switching delay when the latter is varied by 
continuously increasing the ramp's rise (and fall) time from 
1-150 ps. 



considerable precessional motion when the stress is 60 
MPa and virtually no ripple and very little precessional 
motion when the stress is 80 MPa. Thus, by increasing 
the stress slightly from 60 MPa to 80 MPa, the switching 
time can be reduced by a factor of four - from slightly 
above 6 ns to slightly above 1.5 ns. 

The strong dependence of switching delay on stress at 
low stress levels is illustrated in Figure [20] which plots 
switching delay as a function of stress for two different 
rise (and fall) times of 1 ps and 150 ps. Notice that 
switching delay increases rapidly with decreasing stress 
in the interval [60 MPa, 70 MPa] but much less rapidly 
at higher stress levels exceeding 80 MPa. This is purely 
a consequence of the complex out-of-plane dynamics of 
the magnetization vector. This shows that any analysis 
which ignores the out-of-plane dynamics, and tacitly as- 
sumes that the motion of the magnetization vector will be 



always constrained to the plane of the nanomagnet since 
Nd-xx 3> Nd-yy, Nd-zz, will not only be quantitatively 
wrong, but qualitatively wrong as well. 

The cross over between the two curves in Figure [201 
was explained in the context of Terfenol-D and is not 
repeated here. This too is a consequence of the out-of- 
plane dynamics. 



2. Switching delay and energy dissipation 

Figure [5T] shows the energy dissipated in flipping the 
magnetization of the nickel/PZT multiferroic nanomag- 
net as a function of the switching delay. The latter is 
varied by varying the applied stress between 60 and 107 
MPa with fixed rise (and fall) time of 150 ps for the stress 
ramp. The energy dissipated internally in the nanomag- 
net due to Gilbert damping and the 'CV 2 ' energy dissi- 
pated in the external circuit are shown separately. Both 
dissipation components decrease smoothly with increas- 
ing switching delay, implying that the average power dis- 
sipated during switching decreases rapidly with increas- 
ing delay. Both tend to saturate as the switching delay 
becomes longer. 

In Figure EU note that the 'CV 2 ' energy dissipated in 
the switching circuit is 1-2 orders of magnitude higher for 
nickel than for Terfenol-D for the same switching delay. 
Since the voltage V is proportional to stress, the l CV 2 ' 
energy is quadratically proportional to stress. Since the 
magnetostrictive coefficient of Terfenol-D is considerably 
higher than that of nickel, Terfenol-D requires much less 
stress to generate the same stress anisotropy energy and 
hence requires much less stress to switch. This results in 
a significant reduction of 'CV energy dissipation in the 
case of Terfenol-D. 

The total energy dissipation is however dominated by 
the energy Ed dissipated internally in the nanomagnet 
due to Gilbert damping. This energy is actually smaller 
in nickel than in Terfenol-D, because the Gilbert damp- 
ing constant of nickel is more than twice smaller than 
that of Terfenol-D. As long as we are switching adiabat- 
ically, Ed will be the primary source of dissipation, and 
in that case, the material with the lower Gilbert damp- 
ing constant will be superior since it will reduce total 
dissipation. The situation may change completely if we 
are switching abruptly. In that case, the 'CV 2 ' energy 
may very well be the major component of dissipation. If 
that happens, then the material with the larger magne- 
tostrictive coefficient will be better since it will need less 
stress to switch and hence less voltage and less 'CV 2 ' dis- 
sipation. In other words, nickel is better than Terfenol-D 
(from the perspective of energy dissipation) when switch- 
ing is adiabatic, but the opposite may very well be true 
when the switching is abrupt. 

If speed is the primary concern, then what matters 
is the product of the magnetostrictive coefficient and the 
Young's modulus. Since the maximum strain that can be 
generated in the magnetostrictive layer is fixed and deter- 
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mined by the PZT layer, the maximum stress anisotropy 
energy that can be generated in the magnet depends on 
the aforesaid product. The higher the stress anisotropy 
energy is, the faster will be the switching. Since the prod- 
uct is higher for Terfenol-D than for nickel, the Terfenol- 
D/PZT multiferroic switches faster. It can be switched 
in sub-ns by the maximum strain generated in the PZT 
layer, but that same strain cannot switch a nickel/PZT 
multiferroic in less than 1 ns. 

Figure [22] shows the energy dissipated as a function 
of switching delay when the latter is varied by varying 
the ramp's rise (and fall) time between 1 ps and 150 ps, 
while holding the stress constant at 107 MPa. The l CV 2 ' 
component decreases with increasing rise (and fall) time 
because switching becomes increasingly adiabatic. How- 
ever, the dependence of Ed on switching delay is more 
complicated. When the rise and fall of the ramp is fast, 
ripples appear. That causes the energy dissipation to 
increase with increasing delay contrary to expectation. 
This corresponds to the segment of the curve between 
switching delays of 0.9 to 1.9 ns which corresponds to 
rise times between 1 and 120 ps. At slower rise times be- 
tween 120 and 150 ps, the ripples subside and the energy 
dissipation Ed decreases with increasing delay. This cor- 
responds to the segment between 0.9 ns and 1 ns switch- 
ing delay. 



C. Cobalt 

Cobalt has a negative magnetostrictive coefficient that 
is similar to nickel's. Therefore, we will need a tensile 
stress to initiate magnetization rotation away from the 
easy axis. Its Gilbert damping constant is however small- 
est among the three materials considered (see Table QJ 
and hence we expect it to be least dissipative. For the di- 
mensions of the nanomagnet chosen, the minimum stress 
that we will need in a cobalt/PZT multiferroic to switch 
is 56 MPa, while the maximum stress that can be gen- 
erated by the 500 ppm strain in the PZT layer is 104.5 
MPa. 



1. Ramp rate and switching delay 

As before, Equations (fT9l and (|2"0"|) are solved numer- 
ically to find the values of 6(t) and 4>(t) at any given 
instant t for the cobalt/PZT nanomagnet. 

a. Fast ramp: Figure [53] shows the magnetization 
dynamics of a cobalt/PZT multiferroic nanomagnet when 
the stress is ramped up linearly in time from to the max- 
imum possible value of 104.5 MPa in 1 ps. The in-plane 
and out-of-plane dynamics of the magnetization vector 
are very similar to the case of Terfenol-D or nickel, ex- 
cept now we see even more ripples and even more out- 
of-plane precession. Cobalt shows the most prominent 
ripples because it has the smallest Gilbert damping con- 
stant among all three ferromagnets considered. Hence, 



Dynamics of e with time (104.5 MPa stress, 1 ps ramp) 
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(b) 

FIG. 23. Magnetization dynamics in the cobalt/PZT multi- 
ferroic nanomagnet. The stress is ramped up from to 104.5 
MPa in 1 ps: (a) polar angle 9 versus time, and (b) the tra- 
jectory traced out by the tip of the magnetization vector in 
three-dimensional space while switching occurs, i.e. during 
the time 9 changes from 179° to 1°. 



the precessional motion is least damped. 

b. Slow ramp: Figure [24] shows that a slow ramp 
rate (150 ps rise (and fall) time) quenches the preces- 
sion of the magnetization vector. This happens because 
a slower ramp causes less out-of-plane excursion of the 
magnetization vector and hence less precessional motion. 
This reduces the switching delay from 1.6 ns in the case 
of 1 ps rise (and fall) time to 0.85 ns for 150 ps rise (and 
fall) time. 

Figure [21] shows switching delay as a function of the 
ramp's rise (and fall) time for various stresses. At all 
stress levels - not just at high stress levels as in the case of 
nickel - the switching delay decreases with increasing rise 
(and fall) time because of the ripples and the underlying 
precession that are generated when the stress is ramped 
up very fast. Unlike in the case of nickel, this happens 
even at small stress levels in cobalt since cobalt has the 
smallest Gilbert damping constant among the three and 
is hence most susceptible to precession. 
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Dynamics of e with time (104.5 MPa stress, 150 ps ramp) 
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Dynamics of magnetization axis (104.5 MPa stress, 150 ps ramp) 
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FIG. 24. Magnetization dynamics in the cobalt/PZT multi- 
ferroic nanomagnet. The stress is ramped up from to 104.5 
MPa in 150 ps: (a) polar angle 6 versus time, and (b) the 
trajectory traced out by the tip of the magnetization vector 
while switching occurs. 
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Ramp time vs. switching delay (cobalt) 
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FIG. 25. Switching delay of the cobalt/PZT nanomagnet as 
a function of the rise (or fall) time of the ramp, with the 
magnitude of stress as a parameter. 
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FIG. 26. Switching delay versus stress for the cobalt/PZT 
multiferroic nanomagnet for two different ramp rise (and fall) 
times of 1 ps and 150 ps. 



2. Switching delay and energy dissipation 

Figure [55] shows the dependence of switching delay on 
stress for two different ramp rise (and fall) times of 1 
ps and 150 ps. As expected, the switching delay has a 
very weak dependence on the rise (and fall) time at small 
stress levels when not too much ringing occurs. However, 
at high stress levels, the switching delay increases consid- 
erably for the shorter ramp time since that causes pro- 
longed ringing associated with precessional motion. This 
is a more pronounced effect for cobalt than for nickel or 
Terfenol-D since cobalt has the smallest Gilbert damp- 
ing constant among all three that makes it particularly 
susceptible to precessions. 

Figure [57] shows the energy dissipated in flipping the 
magnetization as a function of switching delay when the 
latter is varied by varying the stress between 56 MPa 
and the maximum possible 104.5 MPa. Both the in- 
ternal energy dissipation Ed and the energy dissipated 
in the switching circuit 'CV 2 ' decreases with increas- 
ing delay showing that the average power dissipated de- 



stress: 56-104.5 MPa, Ramp time: 150 ps (cobalt) 
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FIG. 27. Energy dissipated in flipping the magnetization 
of the cobalt/PZT multiferroic nanomagnet as a function of 
switching delay for a ramp rise (and fall) time of 150 ps. This 
range of switching delay corresponds to a stress range of 56 
MPa to 104.5 MPa. The energy dissipated in the nanomag- 
net due to Gilbert damping and the energy dissipated in the 
external switching circuit ('CV 2 ') are shown separately. 



16 



Stress: 104.5 MPa, Ramp Time: 1-150 ps (cobalt) 




1 1.2 14 1.6 
Switching delay (ns) 

FIG. 28. For a fixed stress of 104.5 MPa, energy dissipated 
in flipping the magnetization of the cobalt/PZT nanomagnet 
as a function of switching delay when the latter is varied by 
continuously increasing the ramp's rise time from 1-150 ps. 



creases quite rapidly with increasing delay. Overall, the 
energy dissipation is slightly less here than in nickel be- 
cause cobalt has the lower Gilbert damping constant and 
slightly higher saturation magnetization (see Equation 
([2"5]) . Once again, the dissipated energies saturate at long 
delays. 

Figure [28] plots the energy dissipations as a function of 
the switching delay when the latter is varied by mono- 
tonically varying the rise time of the ramp from 1 to 150 
ps, while holding the stress constant at 104.5 MPa. For 
the Ed plot, there are two segments. In the first segment 
spanning the range of switching delay between 0.75 ns 
and 1.8 ns, the ramp's rise (and fall) time was varied 
from 1 ps to 50 ps. In this interval, the switching delay 
decreases with increasing rise (and fall) time because of 
the ripple effect (see Figure [231 and the 104.5 MPa stress 
curve). In this interval, the energy dissipated internally 
in the magnet Ed goes up with increasing ramp time 
since we are switching faster, but the 'CV 2 ' energy dis- 
sipated in the external circuit goes down since switching 
is becoming more adiabatic (rise time is longer)^. The 
second segment spanning the switching delay range be- 
tween 0.75 and 0.85 ns corresponds to rise (and fall) time 
between 50 ps and 150 ps. In this range of rise (and fall) 
times, the switching delay increases with increasing rise 
(and fall) time because the ripples subside (see Figure [25] 



again and the 104.5 MPa stress curve). In this range, 
Ed goes down with increasing ramp time since we are 
switching slower. The 'CV 2 ' energy dissipation of course 
always decreases with increasing rise (and fall) time since 
the switching becomes increasingly adiabatic. 



V. DISCUSSIONS 

We have analyzed the switching dynamics in a mul- 
tiferroic nanomagnet consisting of a PZT layer and a 
magnetostrictive layer subjected to time-varying stress. 
The stress is ramped up linearly in time with different 
rates or rise (and fall) times. Three different materials 
(Terfenol-D, nickel, cobalt) were considered for the mag- 
netostrictive layer. They show different behavior because 
of different material parameters (Youngs' modulus, mag- 
netostrictive coefficient and Gilbert damping). 

For the type of magnets chosen (materials and dimen- 
sions) , the minimum switching delay that we have found 
is 252 ps obtained with Terfenol-D using a rise (and fall) 
time of 60 ps for a compressive stress of 40 MPa. The 
corresponding 'CV 2 ' energy dissipation in the switching 
circuit is 30 kT and the energy dissipated internally in 
the nanomagnet due to Gilbert damping is Ed — 777 kT. 
For nickel, the minimum switching delay is 930 ps ob- 
tained with a stress of 107 MPa and a rise (and fall) time 
of 120 ps. In this case, the l CV 2 ' energy dissipation in 
the switching circuit is 15 kT and the energy dissipated 
internally in the nanomagnet due to Gilbert damping is 
Ed = 79 kT. Finally, for cobalt, the minimum switching 
delay is 723 ps obtained with a stress of 104.5 MPa and 
a rise (and fall) time of 50 ps. The 'CV 2 ' energy dissi- 
pation in the switching circuit is 36 kT and the energy 
dissipated internally in the nanomagnet due to Gilbert 
damping is Ed — 75 kT. 

All this shows that it is possible to switch multiferroic 
nanomagnets in less than 1 ns while dissipating energies 
of ^100 kT. This range of energy dissipation is far lower 
than what is encountered in spin transfer torque based 
switching of nanomagnets with the same switching de- 
lay. Thus, strain induced switching of magnets has the 
potential to emerge as the preferred method of switch- 
ing magnetic bits in magnet-based non- volatile logic and 
memory. 
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